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ABSTRACT 

Measuring the polarization of the prompt 7-ray emission from GRBs can significantly improve our under- 
standing of both the GRB emission mechanisms, as well as the underlying engine driving the explosion. We 
searched for polarization in the prompt 7-ray emission of GRB 041219a with the SPI instrument on INTE- 
GRAL. Using multiple-detector coincidence events in the 100-350 keV energy band, our analysis yields a 
polarization fraction from this GRB of 98 ± 33%. Statistically, we cannot claim a polarization detection from 
this source. Moreover, different event selection criteria lead to even less significant polarization fractions, e.g. 
lower polarization fractions are obtained when higher energies are included in the analysis. We cannot strongly 
rule out the possibility that the measured modulation is dominated by instrumental systematics. Therefore, 
SPI observations of GRB 041219a do not significantly constrain GRB models. However, this measurement 
demonstrates the capability of SPI to measure polarization, and the techniques developed for this analysis. 
Subject headings: polarization, instrumentation: polarimeters, methods: data analysis, techniques: polarimet- 
ric, gamma rays: bursts, gamma rays: observations 



1. INTRODUCTION 

Despite the extensive work in recent years on GRB after- 
glows, the nature of the central driver that powers the burst 
and the prompt 7-ray emission mechanism remain enigmatic, 
as the physics of the afterglow is insensitive to the nature of 
the progenitor once a relativistic fireball is formed. There has 
been different suggestions for the mechanism that powers the 
GRB central engine. In th e models invok i ng merging neu- 
tron stars and ' collapsars' (Wooslev 1993; Paczvnski 1998; 
MacFadven & Wooslev 1999), hydrodynamically dominated 
outflows (jets) transport the bulk GRB kinetic energy. Al- 
ternatively, Poynting-fiux may be the driver for the tra nsport 
of en ergy to large distances tLvutikov. Pariev & Bla ndford 
2003). Synchrotron radiation has traditionally been the 
favor ed emission m echanis m of the promp t 7-ray emis- 
sion ( Meszaros, Rees & Papathanassiou 1994; Tavani 1996; 
iDermer. Chiang & Bottcherll 19991 iLlovd & PetrosianlhoOOl) . 
though competing Compton upscattering and synchrotron- 
self Compton mod els have been put for ward (Liang 1997; 
Meszaios, Rees & Papathanas sioul Il994t IChiang & Dermerl 
Il999t ISari & Esinll2001l IZhang & Meszaros f|2001l): reviews 
of GR B models can be found in Piran ( 1999) and M eszarosl 
l l2001l) . In terms of polarization mod eling, synchro tron radia- 
tion is natura lly a strong candidate (Coburn & Boggsl 120031: 
iGranotl 120031) . but a portion of the p olarized photon sig- 
nal m ay also be Compton up-scattered ( Eichler & Levinson 
2003). A definite measurement of polarization properties 
from the prompt emission of GRBs will probe their anisotropy 
or magnetic field geometry, and thereby help determine the 
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nature of the central engine and the 7-ray emission mecha- 
nism. 

The first detection of the linear polarization from the 
prompt 7-ray emission of a GRB indicated a very hi gh po- 
larization fraction of 80 ± 20 % dCoburn & Bogg sl 12003). 
For this measurement, Reuven Ra maty High Ene rgy So- 
lar Spectroscopic Imager (RHESSI, iLin et al J 1200 2) data of 
GRB 021206 were used. The measurement demonstrated the 
potential for measuring polarization using Compton-scattered 
events between multiple detectors. Using RHESSI, the same 
method is used later to measure the pol arization fraction of 
two X-class solar flares (Boggs. Coburn & Kalemci 2006). 

The large polari zation fraction obtained by 
Coburn & Boggs (2003) resulted in a series of theoret- 
ical work on 7-ray polarization in GRBs and P o ynting 
dominated flows llLyutikov. Pariev & Blandf ord! 120.231 
Nakar, Piran & Wax manl l2003t lE ichler & Levinson 2003; 
branoll l2003t iDail 12004 lLazzati et alJ 12004 . However, 
independent analyses of the RHESSI data by other groups 
were not ab le to confirm this result at the same leve l of 
significance ( Wigg er et aDl2004t iRutledge & Foxl 12004 . so 
that the degree of polarization for GRB 021206 remains 
uncertain. Clearly, more measurements, using different 
instruments and techniques, are required. Re cently, using the 
BATSE instrument on CGRO, IWillis et all (120051) provided 
evidence for large polarization fractions for two bursts, 
GRB 930131 (n > 35%) and GRB 960924 (TI > 50%), 
without strongly constraining the upper limits. In their work, 
the mass model of BATSE, along with a mass model of the 
earth's atmosphere were used, and the polarization fraction 
was determined by analyzing the angular distribution of pho- 
tons that are scattered throug h the earth's a tmosphere. SPI 
(Spectrometer on INTEGRAL, Vedr enne et alJ2003l) and IBIS 
(Imager on Board the INTEGRAL Satellite. lUbertini etall 
2003) instruments on INTEGRAL can also measure the 
polarization fraction and angle of a source using the coinci- 
dence events between detector pairs ( Lei. Dean & Hills fl997l 
Kale mci et alJ 12004. similar to the method employed by 
ICoburn & Boggsl ( 120031) : iBoggs. Coburn & Kalemci ( 120061) 
with RHESSI. 
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FIG. 1 .— The SWIFT light curves of GRB 041219a in (a) 25-50 keV, (b) 50-100 keV and (c) 100-350 keV bands i Barthelm vet all2004l) . and the SPI light 
curves for the singles in the same bands in (d), (e), and (f) respectively. Times are relative to the peak of the precursor. 



In this letter, we discuss methods to measure polarization 
using one of the instruments on INTEGRAL, SPI, and ap- 
ply these methods to measure the polarization properties of 
GRB 041219a, a bright and a long (about 450 s) GRB which 
was detected with the INTEGRAL Burst Alert System (IBAS) 
and with ISGRI on December 19 at 01:43 UT. The burst 
is in the fully coded field of view of both the ISGRI and 
the SPI, and is ^3° off the X-axis, and 155° in azimuth 
from the Y-axis 8 . The ISGRI coordinates are reported as 
RA = 6.1075° and DEC=+62.8349° with an uncertainty of 
2' (iGotzet all 12004). The brightest part of the burst satu- 
rated the available telemetry of INTEGRAL. The long duration 
and brightness allo wed for multi-wavelength campaigns for 
this GRB (Blake & Bloo mll2004albT: ISoderberg & Frail2004l 
ISonoda et al J 120041) . The infrared counterpart locati on is 
given a s RA=6.1153°, and DEC=62.8426° (Blake & Bloom 
2004b ). The burst was also detected wi th SWIFT-BAT 
(Barthelmy et al. 2004; Fenimore et al. 2004). A comprehen- 
sive spectral and temporal analysis of the burst with SPI, 
SWIFT -BAT, and the RXTE ASM is given in lMcBreen et all 
(2006). The SWIFT-BAT and SPI (singles) light curves are 
shown in Fig. [2 These light curves indicate that the spec- 
trum softens as the burst progresses, and also show a precur- 
sor ^250 s before the main peak. 

2. ANALYSIS 



8 In this work, the azimuthal angles are defined in a plane perpendicular to 
the SPI pointing X-axis, and measured with respect to the SPI Y-axis towards 
the SPI Z-axis. 



2.1. SPI and ^-ray polarization 

SPI is a coded-aperture telescope using an array of 19 
coole d germaniu m detectors for high-resolution spectroscopy 
dVedrenne et alJl2003l) . It works in 20 keV - 8 MeV band, 
and has an energy resolution of ~2 keV below 500 keV. The 
fully coded field of view is 16°, and the angular resolution is 
~3°. At the time of the observation, 17 detectors were active 
due to the failures of Detectors 2 and 17 in orbit. If a photon 
deposits all of its energy into one detector, SPI records this as 
a single event. If a photon interacts through Compton scatter- 
ings with energy deposits in more than one detector, the de- 
tector and channel information for each interaction are saved 
into a multiple event (ME). Even though SPI is not primar- 
ily designed for polarization measurements, these ME data 
are inherently sensitive to polarization as linearly polarized 
gamma-rays preferentially scatter in azimuthal directions per- 
pendi cular to their electric polarization vector (Kal emci et alJ 
2004). 

The two main parameters that determine the sensitivity of 
a multi-detector instrument to gamma-ray polarization are the 
effective area to the multiple-detector scatter events, and the 
average value of the polarimetric modulation factor Q, which 
is the maximum variation i n the azimuthal scattering proba- 
bility for polarized photons (tNovickl 19751 iLei. Dean & Hill's] 
1997). This factor is determined by the scattering cross sec- 
tions, 

(dcr± -dcr\\) 



Q- 



(1) 



(da± + da\\) ' 

where da±, da\\ are the Klein-Nishina differential cross 
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FIG. 2. — (a) The ACS MBW counter rate, (b) Detector dead-time, (c) 
Detector germanium saturation rate. Time is the beginning of pointing. 



sections for Compton scattering perpendicular and parallel to 
the polarization direction, respectively, which are functions of 
the incident photon energy and the Compton scatter angle be- 
tween the incident photon direction and the scattered photon 
direction. For a source count rate of S, and fractional po- 
larization of Il 5 , the expected azimuthal scattering angle (</>) 
distribution is 



^ = ^[l-en s cos2(0-77)]. 

dd) Z7T 



(2) 



Therefore, the "signature" of polarization is a 180°-periodic 
modulation in the distribution of azimuthal scattering angles, 
with a minimum at the polarization angle rj. 

2.2. GRB 041219a SPI data 

A first look at the light curve showed that the SPI data for 
this GRB were affected by the telemetr y saturation p roblems 
which also affected the IBIS data (Got z et al.ll2004l) . Fig.[T] 
shows the SPI light curves (right block) for singles for the sum 
of all detectors. Here we define "singles" - all single detector 
events - as the sum of all SE and PSD events in the Integral 
Science Data Center (ISDC) format. As the flux peaks, a sud- 
den drop occurs in the count rate, which was not observed in 
the SWIFT light curves (left block). 

We inspected several housekeeping parameters to verify 
that the origin of the problem is not something other than 
telemetry saturation. The anti-coincidence system rate, the 
dead-time, and the germanium saturation rate for Detector 
are shown in Fig. |2] Even though there is an increase in 
dead-time, the increase is modest, and cannot account for 
the dropouts in SPI light curves. The germanium saturation 
rate shows no significant deviation from the norm that could 
cause a sudden decrease in the count rate. We then compared 
the "raw" and the "prp" on-board time (processed through a 
standard pipeline at the ISDC), and found that the prp times 
have gaps that are approximately multiples of 0.125 s, indi- 
cating that telemetry packets are missing (S. Schanne, private 
communication). Since the time and duration of the gaps are 
known, an approximate light curve can be reconstructed. The 
100-500 keV (total energy) light curve of ME events, cor- 
rected for effective dead-time due to the missing packets, is 
shown in Fig. [5] Characterizing this effective dead-time is 
important in terms of determining the correct background rate 
for the regions with the packet loss problem. 

2.3. MGEANT simulations 
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FIG. 3. — The observed (black histogram) and the reconstructed (red his- 
togram) light curve of ME events in 100-500 keV band. The gaps are treated 
as dead-time. The vertical solid lines separate Rl, R2, R3 regions (see text). 



To determine the polarization fraction for this GRB, we 
need to compare the measured azimuthal scattering angle 
distribution to the expected distribution for an unpolarized 
and a polarized source from this sky location. The only 
method available for performing this comparison is with de- 
tailed Monte-Carlo simulations. The response to a polarized 
source is characte rized by the polarimetric modulation factor, 
Q, discussed in § 12.11 Since Q is energy-dependent, it will 
depend on the energy spectrum of the source. We therefore 
used simul ation s for two purposes, to obtain the spectral pa- 
rameters (§ I2.4> . and to obtain the modulation factors. 

The simulations are performed using MGEANT 
ISturner et alJ [2000), which is a 7-ray instrument simu- 
lation package developed at NASA/GSFC. The MGEANT 
source code allows several beam geometries and spectra to 
be specified at compile time. A highly detailed SPI mass 
model is used as an input to MGEANT. In order to have 
a complet e response, the mass model of the rest of the 
spacecraft (Ferguson et al. 2003) is also included. This mass 
model is the same as the mass model used to create SPI 
response matrices with MGEANT. More information on 
MGEA NT and the complet e mass model we used can be 
found in lSturner etaD J2003I) . 

2.4. GRB 041219a spectrum 

Detailed and precise determination of the GRB 041219a 
spectral parameters is not necessary for this work as Q is not 
strongly dependent on the exact spectral parameters. There- 
fore, a rough determination of the GRB spectrum is adequate 
for our study. To determine the spectral parameters, we first 
obtained the singles count spectrum from the region with no 
packet loss (Rl in Fig. [3}- We determined the background 
for each detector as follows: We took the data from the first 
1000 s from the beginning of the pointing and obtained a 
spectrum. Next we applied two corrections. We fit the back- 
ground light curve with a first order polynomial to take into 
account a small (a few percent) and gradual increase towards 
the GRB. Second, we found the live-times at the background 
region (B/, ve ,,„ le ) and the GRB region (GRBn ve ,j me ), and mul- 
tiplied the background spectrum with GRB\i vetme IBii vet i me . An 
example spectrum after background subtraction is shown in 
Fig. El 

After the photons from the GRB event were isolated and 
spectra for each detector were obtained, the next step is to re- 
produce these spectra with simulations. To perform the sim- 
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FIG. 4. — The top panel shows the singles light curve of Detector 15 in 
30-490 keV band. No dead-time correction is applied. The solid vertical 
lines indicate the region (Rl) for which the spectrum is extracted. The bot- 
tom panel shows the measured (black histogram), and the simulated (purple) 
spectrum. 



ulations, we first modified the original mass model such that 
the detectors 2 and 17 are not active. We obtained longitude 
and latitude of the spacecraft axes using "spipoinf, and used 
the position of the infrared counterpart for the G RB. We ran 
three simulations with the Band Function JBand et alii 1991 
spectrum using (1) a=l.0, (3=2.4, E hr =\10, (2) a=1.0, (3=2.0, 
E br =170, and (3) a=l.0, (3=2.0, E br =200. 

For simplicity, the mass model uses a single mass for all 
detectors, even though in reality the mass of each detector is 
slightly different. This mass distribution causes the largest 
detec tor-to-detector variations in efficiency JSturner et al.l 
120031) . The simulated spectra were corrected for this effect. 
We also applied a correction for dead-time for each detec- 
tor. Apart from these detector-dependent corrections, there 
are also detector independent corrections regarding the photo- 
peak efficiencies and the mask transmis s ion. These were 
also applied as described in Sturn er et alJ d2003b . We found 
that the spectrum with these set of parameters, a=1.0, (3=2.0, 
Ebr=200, best describes the data in Rl. In Fig.|4] we show 
the actual and the simulated spectrum of Detector 15 as an 
example 9 

We also checked the detector distribution of 30-490 keV 
band total singles counts and compared it to the simulated 
distribution. The result is shown in Fig. [5] The simulation 
(shown with dashed lines) reproduced the actual distribution 
well (within a few %) for detectors not shadowed by the mask. 
For the detectors under the shadow of the mask elements, the 
discrepancy is larger if expressed in percentage. Note that the 
number of events in these detectors are much less than the 
open detectors. Finally we compared the total number of ME 
after cuts (see § [3} with the total number of ME from simu- 
lations. We obtained 534 ME from simulations compared to 
543 events from the data, an agreement within 2%. 

Overall, the Band spectrum with a=1.0, (3=2.0, £V=200. 
represents the GRB in Rl, and the simulation reflects the ac- 
tual detector to detector distribution. 

2.5. Modulation Factor 



9 These s pectral results were obtained before the publication of 
IMcBree n et al. 2006. Even though band function parameters are different, 
the effect of small differences in the energy spectrum is not important for 
polarization measurements. 



Observed counts 
Simulated counts 



5 10 15 

Detector Number 

FIG. 5. — Comparison between the actual (solid histogram) and the sim- 
ulated (dashed histogram) number of counts for each detector for singles, in 
30-490 keV band. 



The next step to measure polarization is to obtain Q (see 
Eq. Q by comparing the azimuthal scattering angle distribu- 
tions of non-polarized and 100% polarized photons. Both 
can be created using MGEANT simulations. Determining 
the azimuthal scattering angle requires finding the direction 
of the photon as it scatters from one detector to the other. SPI 
records the energies and the detectors in a ME. But the di- 
rection of the photon cannot be uniquely determined for all 
events. The conservation of energy and momentum in the 
Compton scattering process place limits on the energies de- 
posited in each detector. Assuming full energy deposition in 
two detectors after a single Compton scattering: 



m,,c 



m e c 
Ei+E 2 



= l-cos(6») 



(3) 



where E\ is the energy deposition in the first detector, and 
Ei is the energy deposition in the second detector, and m e 
is the mass of the electron. One can easily show that for 

relatively small total energies (Ei +E2 < E 2 is always 
greater than E\. As the initial energy increases, the num- 
ber of cases with E2 < E\ increases, and finally at m e c 2 =5 1 1 
keV, there is equal probability for either case. For the spec- 
trum of GRB 041219, most of the photons Compton scattered 
from the low energy deposition detector to the high energy 
deposition. Therefore we tag the direction of every photon as 
originating from the lower energy deposition detector. Even 
though some of the interactions will be tagged incorrectly this 
way, the final results should not be affected significantly due 
to the 180° symmetry of the polarization modulations. 

On the other hand, MGEANT simulations provide more 
information than that of the real data. First, in simulations, 
the interaction positions within the individual detectors are 
known. Second, for any incoming photon energy, the direc- 
tion of the photon is also known. We determined azimuthal 
scattering angle distributions for three cases; (a) using the 
actual interaction positions and directions determined by the 
simulation, (b) using the detector center-center angles (pixel- 
lation) and directions determined by the simulation, (c) using 
the center-center angles and directions determined using en- 
ergy depositions. Cases (a) and (b) can only be calculated 
using the simulations, and (c) represents distribution for the 
actual data. 

We obtained the modulati on factors by following the 
method described in lLei. Dean & Hillsl dl997l) . For the simu- 
lated events with 100% polarized photons (n s =l) the modu- 
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FIG. 6. — Simulated azimuthal scattering angle distribution of 100% po- 
larized photons at 200 keV originating at the GRB position. The center-to- 
center pixellated distribution (black fit) has lower modulation amplitude than 
the distribution obtained using interaction positions within the detectors (red 
fit). 
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FIG. 8.— Simulated modulation factors for GRB 041219a (100-350 keV) 
as a function of polarization angle. Interaction positions within detectors are 
used to obtain the points shown with circles. The triangles, on the other hand, 
uses center-to-center angles. Our best-fit determination of the input spectrum 
during Rl is used. 
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FIG. 7. — The simulated modulation factors as a function of energy for 
three cases; diamonds for the hypothetical case of an on-axis GRB (with the 
same spectra of GRB 041219) with 19 detectors on, using the actual direc- 
tions from simulations (case b), filled circles and crosses are for the known 
position of GRB 041219a with 17 detectors, with case b, and c, respectively 
(see text for more explanation on each case). All three cases pixellated. 



lation factor can be obtained by fitting the azimuthal scatter- 
ing angle distribution with a cos2(</>— 77) function (see Eq.[2}. 
However, before doing this, one needs to take into account the 
"response" of the distribution for non-polarized photons. This 
response is obtained by dividing the non-polarized simulated 
azimuthal scattering angle distribution by its average. We di- 
vided the 100% polarized azimuthal scattering angle distribu- 
tion with this response. For the response we use, the modula- 
tion is on top of an average rate. 

Fig- El shows the azimuthal scattering angle distribution of 
100% polarized photons at 200 keV as an example. The am- 
plitude of the modulation with respect to the average gives 
the modulation factor. The pixellation reduces the modula- 
tion factor around 20% (with respect to non-pixellated modu- 
lation) at 200 keV. 

We ran more simulations with mono-energetic photons at 
different energies, with non-polarized and 100% polarized 
photons, with a randomly chosen polarization angle of 45°. 
Then we histogrammed the azimuthal scattering angles using 
three different methods described earlier. Fig. [7] shows the 
distribution of modulation factor as a function of energy for 
different cases. The ~3° off-axis position of the GRB, and 
the reduced number of detectors did not affect the modula- 



tion factors significantly. More importantly, using the energy 
depositions to determine the directions rather than using the 
actual directions has no effect on the modulation factor. 

Finally, by using simulations with the GRB spectrum de- 
scribed above, we determined the modulation factor for Rl 
in 100-350 keV band. We ran the simulations with differ- 
ent polarization angles. The results are shown in Fig. [H] For 
the pixellated case, Q varies between 18% and 21%. The 
average modulation is 20%. With real interaction positions 
and order, the distribution shows a similar behavior, vary- 
ing between 23% and 26% with an average of 25%. Similar 
to monochromatic tests, the pixellation reduces the average 
modulation factor 20% for the GRB position and spectrum. 

2.6. Search for systematic effects 

Since INTEGRAL is not rotating, systematic effects not 
foreseen by the simulations could alter the azimuthal scat- 
tering angle distribution. Even though we apply correc- 
tions to the simulations, there may be systemati c effects re- 
lated t o ME events that are not discussed in ISturner et alJ 
J2003I) . We therefore analyzed some of the data taken at 
the gro und calibration te sts of SPI (the Bruyeres-Le-Chatel 
dataset. lAttie et alJ2003l) to search for systematic effects that 
could affect polarization measurements. The best calibration 
dataset for our purposes is the case with no mask, on-axis, and 
using I33 Ba, which resulted in strong lines at 276.4, 302.9, and 
356.0 keV. Unfortunately, after the runs, it was discovered that 
the source position was slightly offset. A Hostaform plastic 
device inserted in the center of the plastic anti-coincidence 
scintillator shadows Detector 0. Due to the offset, it also 
casted a shadow on Detectors 2 and 3. We decided to ignore 
all interactions between Detector and the surrounding six 
nearest neighbor detectors for this analysis. We also excluded 
the interactions that involve Detectors 2 and 17, as they are no 
longer operational. We only used the interactions for which 
the total energy gives the line energy. We did not apply a dead- 
time correction, but applied a correction factor to account for 
differing detector masses (see § [5}. Then we ran our stan- 
dard histogram procedures to obtain the azimuthal scattering 
angle distribution. Fig. [9] shows this distribution normalized 
by its average. The variations are in the 1 % level. Therefore, 
excluding interactions in Detector 0, the systematic errors in- 
herent to the detectors are of the order of 1%. 
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FIG. 9. — The azimuthal scattering angle distribution of a calibration run 
with no mask. The histogram is divided by its average to obtain variations 
with respect to the mean. The variations are in the 1% level. 



FIG. 10. — The azimuthal scattering angle distribution of events in Region 
1 , and a cos 2(<f> — r]) fit to the data. The solid line is the average (no polariza- 
tion) and the dashed lines show the maximum and the minimum modulation 
for a 100% polarization fraction. 



2.7. Chance coincidences 

An important factor in polarization experiments using co- 
incidence events is the rate of chance coincidences, events 
in two detectors that occur within the pre-determined coinci- 
dence window. The electronic coincidence window is 350 ns 
for SPI, i.e. two events in different detectors that occur within 
350 ns of each other are recorded as a ME. At the peak of the 
outburst, the total count rate in singles for the detectors that 
are not shadowed is ^400 cts/s. When we eliminate pairs that 
will not obey our selection criteria (see § [3j, the maximum 
singles count rate for each detector is ^150 cts/s. Therefore 
the maximum chance coincidence rate at the peak of the out- 
burst is only 0.008 cts/s per detector pair, which is negligible 
for our measurements. 

3. POLARIZATION MEASUREMENT 

Before analyzing the data, we applied three energy cuts to 
the multiple events: the minimum allowed energy for each de- 
tector in a pair is 26 keV, the minimum allowed total energy 
of a pair is 100 keV, and the maximum allowed total energy 
of a pair is 350 keV. The minimum cuts are necessary to en- 
sure that the events are actual Compton events. They also 
cut a significant portion of small Compton scattering angles 
which contribute less to the modulation factor. The maximum 
cut is required for two reasons. First, due to low count rates 
and low modulation factors, including the very high energy 
par t doe s not improve the measurement. Second, as discussed 
in 32.51 as the total energy increases, the number of events 
with incorrect azimuthal scattering angles increases. To ob- 
tain maximum allowed total energy we considered the sig- 
nal/noise ratio of MEs for different energies, their respective 
modulation factors, and finally the fraction of the incorrectly 
tagged events. And finally, we cut all MEs with total energies 
between 184 keV and 201 keV to remove significant number 
of background photons in the prominent Ge line at 198 keV. 

We then defined two pseudo-detectors (PD) for each detec- 
tor pair (i.e. events that scatter from detector to detector 1 
is PD1, and events that scatter from detector 1 to detector 
is PD2, therefore, it is different than the Pseudo-detector def- 
inition of ISDC). We only used the nearest neighbors. Even 
though it is possible to increase the number of angles by us- 
ing non-neighbor detectors, the number of these events are 
too low to include in the analysis. There are 64 PDs after the 
failures of detector 2 and 17. 

We separated the light curve in three regions. Region 1 is 



from the beginning of the burst to time that the packet-loss 
problems began. Region 2 and Region 3 are determined using 
the source and background rates to maximize the source to 
background ratio. These regions are denoted as Rl, R2, and 
R3 in Fig. [3] The analysis is relatively straightforward for 
Rl. For each PD, the background is determined exactly as 
determined for singles; using the first 1000s of the pointing. 
The rate is again corrected for dead-time and evolution. The 
total-background counts are histogrammed into 6 azimuthal 
scattering angles. The total number of source counts is 545, 
and the total number of background counts is 171. 

The simulated, non-polarized events are corrected for mass 
and dead-time. The dead-time for a pair is calculated using 
both the dead-time due to ACS vetos and the detector elec- 
tronics. We multiplied the number of events in PDs with ™'*'" 2 

where m\ and ;«2 are masses of the detectors that form the PD, 
and m avg is the average mass of all detectors. After these cor- 
rections, we histogrammed the simulated data exactly as we 
histogrammed the real data. To obtain the polari zation frac- 
tion, w e followed the method descri bed i n dLei. Dean & Hills 
1997), and also discussed in Section l2~3l The resultant distri- 
bution and the cos 2(0- rj) fit is shown in Fig. EH The best fit 
modulation amplitude is QH S =21.3 ± 7.6%, corresponding to 
a polarization angle rj = 48.3° ± 3.8°. The \ 2 for this best fit 
is 2.69 for 3 degrees of freedom (DOF). For comparison, the 
X 2 for the best fit assuming no polarization (fiat distribution) 
is 11.00 for 5 DOF. For polarization angles rj ~45°, we cal- 
culate the polarimetric modulation factor as 2=21.2 (Fig.|8|l. 
Correcting the best fit modulation amplitude for this factor 
yields a best fit polarization amplitude of H s = 100 ± 36%, 
providing no upper bound. 

We tried to obtain better constraints by combining Regions 
1, 2 and 3. For Regions 2 and 3, we determined the addi- 
tional dead-time due to the missing packets, and corrected the 
background according to this dead-time. The remaining of 
the analysis is the same as Region 1 . For the combined case 
(R1+R2+R3), the source counts and background counts are 
839 and 389 respectively. Because of the evolution of the 
GRB spectra, the combined spectrum is slightly softer than 
the spectrum of Rl . We determined that a Band function with 
a=1.15, /3=2.4, Ebr=l%0 fits the overall spectrum well. 

The azimuthal scattering angle distribution for the com- 
bined case is shown in Fig. ^2 The fit shown yields a modula- 
tion amplitude of gIL=20.2±6.7% with a minimum at 45.4° 
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FIG. 11. — The azimuthal scattering angle distribution of events in Regions 
1 , 2 and 3 and a cos 2((j> — rf) fit to the data. The solid line is the average 
(no polarization) and the dashed lines show the maximum and the minimum 
modulation for a 100% polarization fraction. 



causing a decrease in the overall modulation. Unfortunately 
the statistics are not good enough to test this hypothesis. 

Our analysis indicates that systematic effects from the two 
inactive SPI detectors, as determined from pre-fiight calibra- 
tion data, should not significantly affect these polarization 
measurements. We do not have knowledge of any further 
systematic effects in orbit that could affect this polarization 
measurement. Analysis of more GRBs may reveal systematic 
effects, with the potential of distinguishing whether the high 
modulation we measured here was a result of high polariza- 
tion fraction, a systematic effect, or just a chance fluctuation. 

The packet loss problem did not play a big role in constrain- 
ing the polarization parameters for this GRB. Without any 
packet loss, there would have been a sensitivity gain of 15%, 
which would not have significantly affected the upper limit 
determination, and may have placed a slightly more stringent 
lower limit. This exercise showed that SPI has a better chance 
of measuring polarization fraction for harder, longer bursts. 



± 5.2°. The modulation factor at this angle is 20.4%, cor- 
responding to U s = 99 ± 33%. The \ 2 for the cos2(<7j-77) 
fit is 4.68 for 3 degrees of freedom (DOF), whereas the \ 2 
for the flat distribution is 15.10 for 5 DOF. Neither of these 
fits represent the data well, as seen in Fig. \^\ and also in- 
ferred from the \ 2 values. Given our measurement uncertain- 
ties, and assuming an unpolarized (flat) distribution, a simple 
Monte-Carlo simulation yields the chance probability of fit- 
ting a modulation of this amplitude as 1.01 %. The best-fit 
polarization yields a lower red uced y 2 over the fit assuming 
no polarization, with an F-test (Bevinaton & Robinson 1992) 
value of 3.34 with 17.3% chance probability (over a flat dis- 
tribution). 

4. DISCUSSION 

We have demonstrated techniques to measure polariza- 
tion of the prompt 7-ray emission of a GRB in the field-of- 
view of SPI on INTEGRAL. However, for GRB 041219a, we 
have not strongly constrained models for the emission mech- 
anism nor the central engine. ME count rate is not high 
enough for statistically significant measurements. For com- 
parison, the RHESSI solar flare polarization measurements 
utilize approximatel y 6500 coun t s and 16000 counts for two 
flares lBoggs. Coburn & Kalemcil d2006). and we use only 839 
source counts for GRB 04 1 2 1 9a. 

Another problem is the dependence of the polarization frac- 
tion on energy cuts. The quoted numbers in this work are for 
the cases with the largest polarization fractions with the high- 
est F-test values compared to a flat distribution. However, 
choosing different energy bands for minimum and maximum 
energies yields lower polarization fractions. For example, us- 
ing 500 keV as the maximum total energy yields a polariza- 
tion fraction of 65 ± 31%, with a minimum at 52°. We obtain 
significantly lower polarization fractions if the 198 keV Ge 
line is not filtered out. If we use energies up to 500 keV and 
do not cut the 198 keV line, the polarization fraction is 55 ± 
30%. The polarization angle may be changing with energy, 



5. CONCLUSION 

We have used data from SPI on INTEGRAL, an instrument 
not intended for polarization studies, and tried to constrain 
polarization parameters of GRB 041219a a long and bright 
GRB in the field of view. The distribution of azimuthal scat- 
tering angles from GRB 041219a is better represented by a 
polarized source compared to a non-polarized source, but with 
low statistical significance. Due to large uncertainties, we 
have not strongly constrained models for the emission mech- 
anism, nor for the central engine. In order to do so, future 
soft 7-ray missions with polarization sensitivity should nec- 
essarily aim for an ability to measure polarizations at the 5- 
10% level, preferably in 3-4 neighboring energy windows, 
and also 2—4 intervals spanning a burst's duration. These re- 
quirements would render a 7-ray polarimeter capable of ex- 
ploring energy-dependent polarization dependence around the 
v — F v peak, and also temporal evolution of both the angle and 
degree of polarization. Knowledge of such source character- 
istics can realistically discriminate between some suggested 
radiation mechanisms and different model geometries that are 
currently being contemplated in the GRB literature. 



E.K. is supported by the European Commission through a 
FP6 Marie Curie International Reintegration Grant (INDAM, 
MIRG-CT-2005-017203). E.K. also acknowledges NASA 
grant NAG5-13142 and partial support of TUBlTAK. M.B. 
acknowledges NASA grant NNG06GB81G for support. We 
thank Steve Sturner and Chris Schrader at GSFC for provid- 
ing the BLC data, and their helps in MGEANT analysis. We 
also thank Volker Beckmann for his helps on the SPI GRB 
analysis. We thank Stephane Schanne at SACLAY for clar- 
ifying telemetry problems, and Pierre Dubath at ISDC for 
clarifying dead-time routines. We acknowledge Mark Kippen 
and Mark McConnell for the "GLEPS" package in MGEANT. 
E.K. thanks Dave Willis and Roland Diehl at MPE for useful 
discussions. 



REFERENCES 

Artie, D., et al., 2003, A&A, 41 1, L71 Blake, C, & Bloom, J. S„ 2004b, GRB Circular Network, 2872 

Band, D., et al., 1993, ApJ, 413, 281 Boggs, S. E., Coburn, W., & Kalemci, E„ 2006, ApJ, 638, 1 129 

Barthelmy, S., et al., 2004, GRB Circular Network, 2874 Chiang, J., & Dermer, C. D., 1999, ApJ, 512, 699 

Bevington, P. R., & Robinson, D. K., 1992, Data reduction and error analysis Coburn, W., & Boggs, S. E., 2003, Nature, 423, 415 

for the physical sciences, New York: McGraw-Hill, Icl992, 2nd ed.) Dai, Z. G., 2004, ApJ, 606, 1000 

Blake, C, & Bloom, J. S., 2004a, GRB Circular Network, 2870 Dermer, C. D., Chiang, J., & Bottcher, M., 1999, ApJ, 513, 656 



8 



Eichler, D., & Levinson, A., 2003, ApJ, 596, L147 
Fenimore, E., et al., 2004, GRB Circular Network, 2906 
Ferguson, C, et al., 2003, A&A, 411, L19 

Gotz, D., Mereghetti, S., Shaw, S., Beck, M., & Borkowski, J., 2004, GRB 

Circular Network, 2866 
Granot, J., 2003, ApJ, 596, L17 

Kalemci, E., Boggs, S. E., Wunderer, C. B., & Jean 2004, in Proceedings 
of the 5th INTEGRAL Workshop, "The INTEGRAL Universe", Munich, 
Germany, Eds. V. Schonfelder, G. Lichti, & C. Winkler, 859 
Lazzati, D., Rossi, E., Ghisellini, G., & Rees, M. J., 2004, MNRAS, 347, LI 
Lei, F., Dean, A. J., & Hills, G. L., 1997, Space Science Reviews, 82, 309 
Liang, E. P., 1997, ApJ, 491, L15+ 
Lin, R. P., et al., 2002, Solar Phys., 210, 3 
Lloyd, N. M., & Petrosian, V., 2000, ApJ, 543, 722 
Lyutikov, M., Pariev, V. I., & Blandford, R. D., 2003, ApJ, 597, 998 
MacFadyen, A. I., & Woosley, S. E., 1999, ApJ, 524, 262 
McBreen, S., Hanlon, L., McGlynn, S., McBreen, B., Foley, S., Preece, P., 

von Kienlin, A., & Williams, O. R., 2006, A&A, 455, 433 
Meszaros, P., Rees, M. J., & Papathanassiou, H., 1994, ApJ, 432, 181 
Meszaros, P., 2001, Science, 291, 79 

Nakar, E., Piran, T, & Waxman, E., 2003, Journal of Cosmology and Astro- 
Particle Physics, 10, 5 



Novick, R., 1975, Space Science Reviews, 18, 389 

Paczynski, B., 1998, ApJ, 494, L45+ 

Piran, T., 1999, Physics Reports, 314, 575 

Rutledge, R. E., & Fox, D. B., 2004, MNRAS, 350, 1288 

Sari, R., & Esin, A. A., 2001, ApJ, 548, 787 

Soderberg, A. M., & Frail, D. A., 2004, GRB Circular Network, 2881 
Sonoda, E., Maeno, S., Matsuo, Y, & Yamauchi, M., 2004, GRB Circular 
Network, 2882 

Sturner, S. J., Seifert, H., Shrader, C, & Teegarden, B. J., 2000, in American 

Institute of Physics Conference Series, 814 
Sturner, S. J., et al., 2003, A&A, 411, L81 
Tavani, M., 1996, Physical Review Letters, 76, 3478 
Ubertini, P., et al., 2003, A&A, 411, L131 
Vedrenne, G., et al., 2003, A&A, 41 1, L63 

Wigger, C, Hajdas, W., Arzner, K., Giidel, M., & Zehnder, A., 2004, ApJ, 
613, 1088 

Willis, D. R., et al., 2005, A&A, 439, 245 

Woosley, S. E., 1993, ApJ, 405, 273 

Zhang, B., & Meszaros, P., 2001, ApJ, 559, 110 



